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Abstract 



o 

D 

Q ■ The two recent density-dependent versions of the finite-range M3Y interaction (CDM3Yn and 

M3Y-Pn) have been probed against the bulk properties of asymmetric nuclear matter (NM) in 
the nonrelativistic Hartree Fock (HF) formalism. The same HF study has also been done with 
the famous Skyrme (SLy4) and Gogny (D1S and DIN) interactions which were well tested in the 
nuclear structure calculations. Our HF results are compared with those given by other many-body 
calculations like the Dirac-Brueckner Hartree-Fock approach or ab-initio variational calculation 



using free nucleon-nucleon interaction, and by both the nonrelativistic and relativistic mean-field 

O ' 

studies using different model parameters. Although the two considered density-dependent versions 



of the M3Y interaction were proven to be quite realistic in the nuclear structure or reaction studies, 
they give two distinct behaviors of the NM symmetry energy at high densities, like the Asy-soft 

ON ' 

and Asy-stiff scenarios found earlier with other mean-field interactions. As a consequence, we 
obtain two different behaviors of the proton fraction in the /^-equilibrium which in turn can imply 

X. 

two drastically different mechanisms for the neutron star cooling. While some preference of the 
Asy-stiff scenario was found based on predictions of the latest microscopic many-body calculations 
or empirical NM pressure and isospin diffusion data deduced from heavy-ion collisions, a consistent 
mean-field description of nuclear structure database is more often given by some Asy-soft type 
interaction like the Gogny or M3Y-Pn ones. Such a dilemma poses an interesting challenge to the 
modern mean-field approaches. 
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I. INTRODUCTION 



The determination of the nuclear equation of state (EOS), which is of vital importance 
for the nuclear astrophysics, has been a central object of numerous studies of heavy ion 
(HI) collisions for the last two decades. If the efforts in the early 90 's were concentrated 
on the determination of the incompressibility K of symmetric nuclear matter (NM), for 
different types of the EOS are usually distinguished by different K values recent studies 
of nuclear reactions involving unstable nuclei lying close to the neutron or proton driplines 
provide us with a unique opportunity to learn about the EOS of asymmetric NM which has 
a large difference between the neutron and proton densities Qj. The knowledge about the 
NM symmetry energy (a key ingredient in the EOS of asymmetric NM) is vital not only in 
studying the dynamics of HI collisions involving radioactive nuclei and/or their structure, 



but a 



,3- 



so in studying the neutron star formation or the r-process of stellar nucleosynthesis 



4], |j] . In particular, the most efficient process of the neutron star cooling, the so-cal 



0,3, 



cd 



direct Urea process in which nucleons undergo direct beta (and inverse-beta) decays 
can take place only if the proton-to-neutron ratio exceeds 1/8 or the proton fraction x > 1/9 
in the /3-equilibrium. The latter is entirely determined from the density dependence of the 
NM symmetry energy S(p) by the following balance equation 

hc(37r 2 px) 1/3 = 4S(p)(l - 2x), (1) 

where p is the total NM density. It is of fundamental importance whether the direct Urea 
process is possible or not. If the x value cannot reach the threshold for the direct Urea 
process, then the neutron star cooling should proceed via the indirect or modified Urea 
process which has a reaction rate of 10 4 ~ 10 5 times smaller than that of the direct Urea 
process and implies, therefore, a much longer duration of the cooling process 6|,LZl]. Although 
a recent test 9J of the microscopic EOS against the measured neutron star masses and flow 
data of HI collisions has shown that the direct Urea process is possible in some cases if the 
predicted neutron star mass is above a lower limit of 1.35 ~ 1.5 solar mass (M ), the overall 
cooling time of a neutron star is still unknown as yet 5J due to the uncertainty about the 
high-density behavior of S(p). We also note here a very strong influence by the neutron 
superfluidity in the inner crust of neutron star 0, [iq|, U], where the superfluid effects can 



reduce the cooling time by a factor of 3 ~ 4. In the same direction, the cooling time might 
be further shortened by the possible phase transitions to quark matter or pion condensate 



occurring in the core of neutron star and leading, therefore, to a much higher central density 
and a smaller star radius Q]. 

Microscopic studies of the EOS of asymmetric NM have been performed in both nonrela- 
tivistic and relativistic nuclear many-body theories, using realistic two-body and three-body 
nucleon-nucleon (NN) forces or interaction Lagrangians (see more details in recent reviews 
0, Q])- These many-body studies have shown the important role played by the Pauli block- 
ing as well as higher-order NN correlations in the G-matrix used to generate the NM binding 
energy at different densities. These medium effects are considered as physics origin of the 
density dependence introduced into various effective NN interactions used presently in the 
nonrelativistic mean-field approaches. Among different versions of the effective NN interac- 
tion, very popular choice is the so-called M3Y interaction which was originally constructed 
by the Michigan State University (MSU) group to reproduce the G-matrix elements of the 



Reid 



13| and Paris [14j NN potentials in an oscillator basis. The use of the original density 



independent M3Y interaction in the Hartree-Fock (HF) calculation of nuclear matter 15] 
has failed to saturate NM, leading to a collapse at high densities. Since the HF method 
is the first order of many-body calculation, some realistic density dependences have been 



introduced into the original M3Y interaction 



16] to effectively account for the higher-order 



NN correlations which cause the NM saturation. During the last decade, different density 



dependent versions of the M3Y interaction 
ric and as ym metric NM [l?], Q, 



19 



20 



rave been used in the HF calculations of symmet- 



21 



22j,l23 



24j , in the mean- field studies of nuclear 



structure [22J, 123|, 124j as well as in numerous folding model studies of the nucleon-nucleus 
and nucleus-nucleus scattering 17, isl, Ji], 20, [25]. In view of these studies, it is now highly 
desirable to have a realistic version of the effective NN interaction for consistent use in the 
mean-field studies of NM and finite nuclei as well as in the nuclear reaction calculations. As 
an exploratory step, we perform in the present work a systematic HF study of NM using 
two density- dependent versions of the finite-range M3Y interaction, the CDM3Yn interac- 



tions which have been used main 
nucleus-nucleus scattering 
parametrized by Nakada 



20 



22, 3, 24] 



y 



in the folding model studies of the nucleon-nucleus and 



25| and M3Y-Pn interactions which have been carefully 



241 ] for use in the mean-field studies of nuclear structure. 



For comparison, the same HF study is also performed with some realistic versions of the 



Gogny [26j, [27j and Skyrme 



2n 



interactions. 
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II. HARTREE FOCK CALCULATIONS OF NUCLEAR MATTER 



Like other mean-field approaches, we consider in the present HF study a homogeneous 
spin-saturated NM at zero temperature which is characterized by given values of neutron and 
proton densities, p n and p p , or equivalently by its total density p = p n + p p and its neutron- 
proton asymmetry 5 = (p n — p p )/(p n + Pp)- With the direct (i>d) and exchange (t>Ex) parts 
of the interaction determined from the singlet- and triplet-even (and odd) components of 
the central NN force, the total NM binding energy is determined within the HF formalism 
as 

E = E^ n H — [< kar, k' o'r^v^kaT, k'a'r' > + < kar, k' o't'\v-exW or, ka'r' >], 



kar k'a'r' 



(2) 



where | kar > are the ordinary plane waves. The nuclear matter EOS is normally classified 
by the NM binding energy per particle which can be expressed as 



j(p,6) = j(p,5 = 0) + S(p)5 2 + O(5*) + ... 
The NM pressure P and incompressibility K are then calculated as 



(3) 



P(p, 6) = p' 



d_ 

dp 



K(p,5) = 9p< 



dp 2 



(4) 



The contribution of 0(S A ) and higher-order terms in Eq. ([3]), i.e., the deviation from the 



parabolic law was proven to be negligible 



17 



291 ] and the most important physics quantity 



is, therefore, the NM symmetry energy S(p) which is the energy required per nucleon to 



29| | . The value of S(p) at 



change the symmetric NM into a pure neutron matter 
the saturation density of symmetric NM, p pa 0.17 fm -3 , is also known as the symmetry 
energy coefficient J = S(pq) which has been predicted by numerous many-body calculations 



to be around 30 MeV 
around p [33[ as 



17 



29 



30 



3l|. The NM symmetry energy S(p) is often expanded 



S(p) = J 



Lfp 

3 



Po 



K sjm i p — po 
Po 



(5) 



where L and K sym are the slope and curvature parameters at p . As discussed above, the 
knowledge about the density dependence of S(p) is of vital importance in studying the 
neutron star formation and has been, therefore, a longstanding goal of many NM studies 



4 



using either microscopic or phenomenological models. The main method to probe S(p) 
associated with a given (mean-field) interaction is to test this interaction in the simulation 
of HI collisions using transport and/or statistical mode 



the structure studies of neutron-rich nuclei 



24 



27 



32 



3 



34. 



35, 



11 



3t, 



37. 



43 



38 



39] or in 



44J. Based on 



the constraints set by these studies using the latest experimental data, some extrapolation 
is then made to draw conclusions on the low- and high-density behavior of S(p). We show 
below that such conclusions still remain quite divergent in some cases. 

In this work we study asymmetric NM at different neutron-proton asymmetries S us- 



ing two different sets o 
CDM3Yn (n = 3,4,6) 



re density-dependent M3Y interaction named, respectively, as 
and M3Y-Pn {n = 3,4,5) 24j. Concerning the first set, the 



isoscalar density dependence of the CDM3Y3, CDM3Y4 and CDM3Y6 interactions has 



been parametrized [I8j] to properly reproduce the saturation point of symmetric NM and 
to give K = 217, 228 and 252 MeV, respectively, in the HF approach ©-flU). These in- 
teractions, especially the CDM3Y6 version, have been widely tested in numerous folding 
model analyses of refractive a-nucleus and nucleus-nucleus scattering (see a recent review 
in Ref. 45]). In the present work, the isovector density dependence of the CDM3Yn inter- 



actions is parametrized using the procedure developed in Ref. [20|, so as to reproduce the 
Brueckner-Hartree-Fock (BHF) nuclear matter results for the energy- and density-dependent 
nucleon optical potential (OP) of Jeukenne, Lejeune and Mahaux (JLM) 46]. This isovec- 
tor density dependence is then scaled by a factor ~ 1.3 deduced recently from the folding 
model analysis [jij, Q of the (p, n) reactions leading to isobaric analog states (IAS) in var- 
ious targets. In contrast to the CDM3Yn set, the M3Y-Pn interactions have been carefully 
parametrized by Nakada [24] in terms of the finite-range M3Y interaction supplemented with 
a zero-range density-dependent force, to consistently reproduce the NM saturation proper- 
ties and ground-state (g.s.) bulk properties of double-closed shell nuclei as well as unstable 
nuclei close to the neutron dripline. These latest versions of the M3Y-Pn interaction have 
not been used in the HF study of asymmetric NM and it is, therefore, of interest to probe 
them here. For completeness, the HF calculation (J2J)- (JHI) has also been done with two other 
popular choices of the effective NN interaction: the Gogny (D1S, DIN) 



26. 



271 ] and Skyrme 



SLy4) 28| forces. Among them, the latest version of the Gogny force DIN has been shown 



27] to reproduce the neutron matter EOS better than the older D1S version while still 



giving a good description of the bulk properties of finite nuclei. 
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TABLE I: HF results for the energy En = E/A(pn,5 = 0) and incompressibility K of symmetric 
NM, symmetry energy coefficient J, slope L and curvature K sym parameters of S(p) evaluated at 
p = po using the CDM3Yn, M3Y-Pn, Skyrme (SLy4) and Gogny (PIS, DIN) interactions. Similar 
results give n by the DBHF calculation using the Bonn A interaction [43] and by other mean-field 



studies 



3 



41 



42 



481 ] are also presented for comparison. K T = K sym — 6L. 
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The properties of symmetric and asymmetric NM described in the HF approximation 
using CDM3Yn and M3Y-Pn interactions are summarized in Table [H Since the isovector 
component of the interaction does not contribute to the total energy of symmetric NM, the 
interactions having about the same values of nuclear incompressibility K are expected to 
give similar EOS up to moderate values of the NM density. One can see from the upper panel 
of Fig. [1] that the HF results for the energy of symmetric NM obtained with the considered 
interactions, which give K « 200 ~ 250 MeV (see Tabled]), are quite close to each other 
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FIG. 1: (Color online) Energy per particle E/A of symmetric NM and pure neutron matter calcu- 
lated in the HF approximation d2J)— d3J) using the effective NN interactions given in Tabled Circles 



491 ] by Akmal, Pandharipande and 



are microscopic results of the ab-initio variational calculation 
Ravenhall (APR). 

at densities up to around 2p - At densities above 2p the NM energies calculated with 
the M3Y-Pn interactions are significantly larger than those obtained with the CDM3Yn 
interactions, even though the predicted incompressibilities K are rather close. The main 



difference here is that the M3Y-Pn interactions have been carefully parametrized [24j not 
only to reproduce the saturation properties of symmetric NM like the parameter choice for 
the CDM3Yn interactions [ljsj], but also to give good description of the g.s. shell structure of 
the magic nuclei. For a comparison, we have also plotted in Fig.[T]the microscopic prediction 
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by Akmal, Pandharipande and Ravenhall (APR) 49j] for both symmetric NM and pure 
neutron matter based on the variational chain summation method, using Al8+5w+UIX* 
version of the Argonne NN interaction. 

For the symmetric NM, the HF results given by all considered interactions agree well with 
the APR predictions at NM densities up to 2p , but at higher densities the APR results 
seem to be closer to the CDM3Y6 and DIN curves. For the neutron matter the picture is 
quite different, with the HF results given by the CDM3Yn and Skyrme (SLy4) interactions 
agreeing with the APR results at the NM densities up to 4p , while the neutron matter 
energies given by the M3Y-P??, and Gogny (D1S, DIN) interactions are by a factor of 2 ~ 3 
smaller than the APR results in the same density range (see lower panel of Fig. [T]). It is 
obvious that such a large difference seen in the HF results for the neutron matter energy 
is due to the difference in the isovector parts of the considered interactions. Since the 
isovector density dependence of the CDM3Yn interactions has been parametrized 20| to 
reproduce simultaneousl y th e BHF results for the isospin- and density dependent nucleon 
OP by the JLM group 46| and charge exchange (p, n) data for the IAS excitations, the 
high-density behavior of the neutron matter energy given by the CDM3Yn interactions 
should approximate that given by a BHF calculation of the neutron matter. In this sense, 
a similarity between the HF results given by the CDM3Yn interactions and microscopic 



APR results 



49 | is not unexpected. In contrast to the CDM3Yn interactions, the isovector 



density dependence of the M3Y-Pn, D1S and DIN interactions were carefully fine tuned 
against the structure data observed for the neutron (and proton-) dripline nuclei and it is 
also natural to expect that the EOS of the neutron matter predicted by these interactions 
should be quite realistic. 

The NM pressure @ is straightforwardly evaluated from the NM energy. The HF results 
for the NM pressure given by the present mean-field interactions are compared in Fig. [2] with 
the empirical constraints deduced from the analysis of the collective flow data measured 
in relativistic HI collisions 34|. For the symmetric NM, all considered interactions give 
consistently the NM pressure well within the borders of empirical data at densities up to 
around 4p (see lower panel of Fig. [2]). For the pure neutron matter, the HF results given 
by the CDM3Yn and SLy4 interactions agree nicely with the data, while those given by the 
M3Y-Pn and Gogny interactions are significantly below the data (see upper panel of Fig. [2]). 
Among the two Gogny forces, the D1S version gives a too low pressure in the neutron matter 
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FIG. 2: (Color online) Pressure of pure neutron matter (upper panel) and symmetric NM (lower 
panel) calculated in the HF approximation ©-([I]) using the effective NN interactions given in 



Table HI The shaded areas are the empirical constraints deduced from the HI flow data 



which fails badly in the comparison with the data. This result confirms again that the D1S 
interaction is not suitable for the study of asymmetric NM as found in the previous NM 



studies 



49 



50 



51] . Looking at Fig. 1 of Ref. 27| one might expect that a Gogny-type 



interaction giving a neutron matter EOS steeper than that given by the DIN interaction 
and closer to the Friedman-Pandharipande's curve 50j would improve the HF description 
of the neutron matter pressure. From Eqs. (jBJ) and (jlj) one finds easily that the difference 
observed in the upper part of Fig. [2] is directly related to different density dependences of 
the NM symmetry energy S(p), which in turn is determined by the isovector component of 
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the considered mean-field interaction. 

To show further the difference caused by the isovector component of the interaction, we 
have plotted in the upper panel of Fig. [3] the HF results for the NM symmetry energy S(p) 
given by the considered interactions in comparison with both the empirical data at low 
densities and microscopic APR results. At the saturation density po of symmetric NM all 
the models predict the symmetry coefficient J « 29 ± 3 MeV which agrees reasonably well 



with the empirical values deduced recent 
(p, n) reaction exciting the IAS states [19 



y from the CC analysis of the charge-exchange 



20[ and structure study of the neutron skin in 



medium and heavy nuclei 52j, |53J. In the low-density region (p 0.3 ~ 0.6 pn 



35 



there exist 
and 



37 



some empirical points extracted from the HI fragmentation data analysis 
the S(p) values given by the CDM3Yn interactions are in a very good agreement with these 
data. The low-density S(p) values given by other effective interactions are slightly larger 
than the data but rather close to the microscopic APR result 49[ . Although the studies of 
HI fragmentation data and/or neutron skin thickness have put some constraints on the NM 
symmetry energy S(p) at p < po, its behavior at higher NM densities still remains uncertain. 
In contrast to the CDM3Y?i, SLy4 and APR predictions, the NM symmetry energies given 
by the remaining mean-field interactions reach their maximal values at NM densities around 

see upper 



37 



54, 



55] 



1.5po an d smoothly decrease to negative values at densities approaching 4 
panel of Fig. [3]). These two different behaviors have been observed earlier | 
and are often discussed in the literature as the Asy-stiff (with symmetry energy steadily 
increasing with density) and Asy-soft (with symmetry energy reaching saturation and then 
decreasing to negative values) behaviors. While Sly4 and other Slyn versions of the Skyrme 



28 



551 ] , the Gogny and M3Y-Pn 



interaction could also be assigned to be the Asy-stiff type 
interactions belong definitely to the Asy-soft type. 

If we compare these two behaviors in terms of the NM pressure shown in Fig. [21 then the 
Gogny and M3Y-Pn interactions clearly fail to reproduce the empirical pressure P(p) of the 
pure neutron matter deduced from the HI flow data 34| . While this comparison might allow 
us to conclude that the (Asy-stiff) CDM3Yn interactions have a more appropriate isovector 
density dependence compared with the (Asy-soft) Gogny and M3Y-Pn interactions, it is not 
possible to do so based on the nuclear structure results given by the CDM3Yn interactions. 
Namely, we have performed spherical HF calculations for finite nuclei and found that the 
CDM3Yn interactions give a much worse description of the g.s. properties of light- and 
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FIG. 3: (Color online) Upper panel: NM symmetry energies S(p) calculated with the interactions of 



531 ] , HI fragmentation 



35 



37 



38] 



Table |H Empirical data are taken from the studies of neutron skin 
and (p, n) excitation of IAS states [19I.I20I]; Lower panel: Proton fractions ( 
interactions, in comparison with the threshold xdu f° r the direct Urea process s| . The curves have 
the same notations as in Figs. Q] and [2j 



medium mass nuclei with neutron excess compared with the Gogny or M3Y-Pn interactions. 
A widely adopted procedure so far is to assume that the NM properties predicted by an 
effective NN interaction should be quite realistic if this interaction gives systematically good 
description of different structure properties of finite nuclei, especially, the unstable nuclei 
near the neutron dripline. In this sense, the failure of the Gogny and M3Y-Pn interactions 
in reproducing the empirical NM pressure as well as the unsuccessful use of the CDM3Yn 



11 



interactions in the nuclear structure calculation pose a serious dilemma, which makes it 
difficult to conclude unambiguously about the high-density behavior of the NM symmetry 
energy S(p) from the present HF results. Although there is no accurate systematics available 
for all existing effective NN interactions on their consistent use in the mean-field studies 



of both NM and finite nuclei, an extensive mean- field study of NM by Stone el al. [55] 
using different Skyrme interactions has shown that out of 87 considered versions of the 
Skyrme interaction only 27 versions give the Asy-stiff behavior of the NM symmetry, and all 
the remaining 60 versions give the Asy-soft behavior. It remains, therefore, an interesting 
question whether the fitting procedure to determine parameters of an effective NN interaction 
(from an optimal description of the g.s. properties of finite nuclei and saturation properties 
of symmetric NM), more likely "bends" the NM symmetry energy S(p) curve down to some 
Asy-soft type rather than "raises" it to some Asy-stiff type. We also note in this connection a 



recent systematics by Klupfel el al. [56j on different phenomenological choices of parameters 
in the Skyrme-Hartree-Fock model for a self-consistent description of nuclear structure and 
NM properties, which shows that the range of NM properties still remains quite broad 
despite a large sample of the nuclear ground-state properties used in the parameter fit. 
It is, therefore, clear that the ability of nuclear structure data to constrain the EOS for 
asymmetric NM on the mean-field level is still limited. 

The difference found between the Asy-soft and Asy-stiff scenarios becomes more drastic 
in terms of the proton fractions (CQ) shown in the lower panel of Fig. [3l If the proton fraction 
x = Pp/(Pp + Pn) exceeds a critical value xdu, then the direct Urea (DU) process 

n — > p + e~ + u e , p + e~ — * n + v e (6) 

becomes possible [f], 0]. Since the neutrino and antineutrino momenta are negligible com- 
pared with those of protons, neutrons and electrons, the DU threshold xdu can be estimated 
from the momentum conservation and charge neutrality condition for p, n and e~ only jf], • 
We have plotted in the lower panel of Fig. [3] the averaged DU threshold as a function of 
the NM density taken from Ref. j^j. At low densities xdu ~ 1/9 as found by Lattimer et 
al. js} in the muon-free approximation. At densities above po, the charge neutrality is cor- 
rected by the muon presence which slightly enhances xuv 2| . From the HF results obtained 
with the considered effective interactions, only the proton fractions given by the CDM3Yn 
interactions can reach the DU threshold at moderate densities p ~ 0.6 fm -3 . 
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According to the microscopic APR study 49J, such a central density is reachable in a 
neutron star having mass M ps 1.6M which is well above a lower limit of 1.35 ~ 1.5M 
for the DU process established in Ref. The NM density p ~ 0.6 fm -3 happens also 
to be within the range of average central densities of the neutron star estimated from a 
nuclear EOS with K ps 240 MeV 6| which is quite close to K values given by the CDM3Yn 
interactions (see Tabled]). As a result, the direct Urea process ()6]) is possible if one chooses the 
CDM3Yn interactions for the in-medium NN interaction in the neutron star matter. Such 
a scenario for the DU process is also favored by the HF results for the NM pressure shown 
in Fig. [2] where the CDM3Yn interactions consistently give the best description of empirical 
data for both the symmetric NM and pure neutron matter. In contrast to the CDM3Yn 
interactions, the choice of the (Asy-soft) Gogny or M3Y-Pn interactions would definitely 
exclude the possibility of the DU process because the corresponding proton fractions can 
never reach the DU threshold as shown in Fig. [3j The microscopic APR results obtained with 
the Al8+5w+UIX* version of the Argone NN interaction approach the muon-free threshold 
Xdu ~ 1/9 only at p ~ 0.8 fm -3 . Such a central density can exist only if the neutron star 
mass M ^ 2M (see Fig. 11 of Ref. 49[) and the DU process is, therefore, very unlikely with 
the EOS given by the APR model. In the case of SLy4 interaction, the numerical integration 
of the Tolman-Oppenheimer-Volkov equation in Ref. 55j has shown that the central density 
reached in a neutron star having mass M ps 1.4M is only p ps 0.55 fm -3 . Since the proton 
fraction x given by the Sly 4 interaction is reaching the DU threshold at much a too high 
density of p ps 1.4 fm -3 55j, the DU process is also not possible with the Sly4 interaction. 



The DU process has also been considered in the fully microscopic many-body studies of 
the EOS using realistic free NN interactions and we found it complementary to com- 

pare t he p resent HF results with those of a recent Dirac-Brueckner Hartree-Fock (DBHF) 
study 47J using an improved treatment of the Bonn-A interaction. It can be seen from the 
upper panel of Fig. |4] that the NM symmetry energy curve given by this DBHF study is 
somewhat stiffer than that given by the CDM3Yn interactions. As a result, the proton frac- 
tion estimated from the DBHF results is reaching the DU threshold already at NM densities 
p ps 0.45 fm -3 . Such a critical density for the DU process is higher than that (p ~ 0.37 
fm -3 ) given by the earlier DBHF results (see Fig. 2 of Ref. and it should correspond to a 
star mass above the lower limit of 1.35 ~ 1.5M for the DU process. It is interesting to note 
that the inclusion of three-body forces into the many-body BHF calculations 58( not only 
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essentially improves the description of the saturation properties of symmetric NM but also 



gives a much stiff er NM symmetry energy at high densities (see Fig. 4 of Ref. [58|). in the 
opposite direction from the Asy-soft type interactions. Given highly accurate parametriza- 
tions of the bare NN interaction, these microscopic many-body calculations are practically 
parameter-free and it is natural to assume an Asy-stiff behavior of S(p) which allows both 
the direct and indirect Urea processes to take place during the neutron star cooling. It also 
is highly desirable that results of such a microscopic many-body study can be accurately 
reproduced at the mean-field level using some effective (in-medium) NN interaction which is 
also amenable to the nuclear structure and/or reaction calculations. However, such "micro- 
scopic" mean-field interactions remain technically complicated to construct and most of the 
structure and reaction studies still use different kinds of the effective NN interaction with 
parameters adjusted to the optimal description of structure and/or reaction data. 

An effective NN interaction can be either fully phenomenological like the Skyrme forces 
or partially based on a microscopic many-body approach like the CDM3Yn interactions 
considered above. An interesting alternative approach has been suggested recently by the 
Tiibingen group 48( which considers only the low-momentum (below a cut-off A = 2 fm -1 ) 
part of the bare NN interaction. While this "low k" interaction Vi ow k still describes well 
the NN scattering data up to the pion threshold, the short-range correlations originated 
from the high-momentum components are treated phenomenologically at the mean-field 
level. Namely, the Vi ow k has been supplemented by an empirical density-dependent contact 
(CT) interaction adjusted to reproduce the saturation properties of symmetric NM and the 
empirical symmetry energy J within the HF approximation 48J. This Vi ow k+CT interaction 
was shown to give also a reasonable description of the g.s. properties of some finite nuclei 
including 208 Pb. The NM symmetry energy and proton fraction predicted by the Vi ow k+CT 
interaction are shown in Fig. [4] and they are quite close to those predicted by the CDM3Yn 
interactions shown in Fig. [3j Like the CDM3Yn interactions, the Vi ow k+CT interaction 
should also belong to the Asy-stiff type and allow both the direct and indirect Urea processes 
during the neutron star cooling. Another famous choice of the effective NN interaction is 
;he Skyrme-type momentum dependent interaction (MDI) which has been first parametrized 



391 ] for the transport model simulation of HI collisions. By varying the x parameter of the 



MDI interaction, the experimental data from NSCL-MSU on the isospin diffusion have been 
shown to favor the MDI (x=-l) version which gives the NM symmetry energy nearly linear in 
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FIG. 4: (Color on 
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Upp er panel: NM symmetry energies S(p) given by different many-body 



47], 



48 



491 ]: Lower panel: Proton fraction (p} given by these symmetry 
energies in comparison with the threshold xdu f° r the direct Urea process j^. See more details in 
text. 



the NM density (see Fig.H]or Fig. 1 of Ref. 39]). One can see in Fig.H]that the NM symmetry 
energy S(p) given by the MDI (x=-l) interaction is somewhat stiffer than that predicted by 
the DBHF calculation using the Bonn A interaction 47j. The proton fraction given by this 
(Asy-stiff) MDI (x=-l) interaction is reaching the DU threshold at NM densities p ~ 0.3 
fm -3 and is well above xdu at the typical central density p « 0.5 ~ 0.7 fm -3 of neutron 
star. Since the star mass corresponding to such a central NM density should be larger than 
the lower limit of 1.35 ~ 1.5M for the DU process 9J, the DU process must be possible in 
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this case. The MDI interaction has also been used to describe neutron skin in finite nuclei in 
the Skyrme HF model 



591 ]. and the MDI interaction with x between and -1 was shown to 
reproduce reasonably well the empirical neutron-skin data for 124 ' 132 Sn and 208 Pb. However, 
the situation with the MDI interaction becomes somewhat confused after the new FOPI 
data on the 7t~/tt + ratio measured in central HI collisions at SIS/GSI energies have been 
shown to clearly favor the MDI (x=l) interaction 36j. In terms of symmetry energy, the 
MDI (x=l) interaction belongs to the Asy-soft type (see Fig. H|) like the Gogny or M3Y-Pn 
interactions considered above and it excludes, therefore, the DU process during the neutron 
star cooling. Given experimental evidences favoring both the Asy-stiff and Asy-soft versions 
of the MDI interaction, the behavior of the NM symmetry energy at high densities as well 
as the possibility of the DU process still remain an open question. 

The NM symmetry energy has also been the subject of various relativistic mean field 
(RMF) studies. In the present work, we compare our nonrelativistic HF results with those 

y chosen parameters for the energy-density func- 



40| | has been shown to consistently reproduce the 



of some recent RMF studies using careful 
tional [4y, |4l|, |42j. The G2 parameter set 
g.s. structure of finite nuclei and bulk properties of NM. In particular, the RMF calculations 
using the G2 parameters reproduce very well the empirical pressure for both the symmetric 
NM and pure neutron matter (^(J. Quite interesting is the FSUgold parameter set developed 



by Piekarewicz et al. 



411 ] which has been used to study not only the NM properties and 



g.s. structure of finite nuclei but also the excitation of the nuclear giant monopole resonance 
(GMR). While the FSUgold parameters give a good description of the GMR in 208 Pb, the 
observed GMR excitation energies in Sn isotopes could not be reproduced using this param- 
eter set. In order to improve the RMF description of the NM saturation properties as well as 
the monopole strength distribution in Sn isotopes, a hybrid model of the RMF parameters 



has been developed [42| based on the earlier NL3 model. However, the Hybrid parameters 
turned out to give a worse description of the GMR in 208 Pb compared with that given by the 
FSUgold model and it remains, therefore, difficult to choose between these two parameter 
sets. The RMF results using these parameters are shown in Fig. 0] and one can see that the 
stiffness of the NM symmetry energy is gradually increasing as one goes from FSUgold and 
G2 to the Hybrid results. It has been found by Steiner 57j that the RMF models typically 
have a large symmetry energy and a large proton fraction, and the DU process becomes 
possible at rather low NM densities. This effect can be clearly seen in the Hybrid and G2 



1(3 



results shown in Fig. [4] where the corresponding proton fractions reach DU threshold at the 
NM densities of p ~ 0.24 and 0.32 fm -3 , respectively. We note further that S(p) predicted 
in the Hybrid model is very close to the RMF result by Klahn et al. 9[ using the NLper 
parametrization where the proton fraction is reaching the DU threshold at p ~ 0.28 fm -3 . 
The behavior of the proton fraction predicted by the FSUgold model is somewhat different 
from those predicted by the Hybrid and G2 models. Namely, it approaches the muon-free 
threshold xbv ~ 0.11 only at p m 0.8 fm -3 , like the microscopic APR result. With the 
predicted maximum neutron star mass of M ~ 1.72M , the FSUgold model has been shown 



in Ref. 41| to allow partially the DU process in the neutron star cooling. However, if we 
adopt the averaged DU threshold taken from Ref. Q which takes into account the muon 
presence at high densities, then the DU process is unlikely in this case because the proton 
fraction predicted by the FSUgold model seems to saturate at densities p ^ 0.8 fm -3 , around 



a value of x ~ 0.11 (see lower panel of Fig. H]) like the APR results [49| discussed above. 

Although the experimental evidences are still divergent with respect to the Asy-stiff 
and Asy-soft type mean-field interactions, it is of interest to further explore the difference 
between these two groups in terms of the NM incompressibility. Using the general definition 
©-(JI]), the NM incompressibility K(p) can be written explicitly in terms of the isoscalar 
{Kq) and isovector (Ki) parts as 

K(p, 5) = K (p) + Ki(p)5 2 + 0(5 4 ) + ... (7) 



where K (p) = 9p 







2 



dp 2 



K M ^ 9 -^. (8) 



It is clear from Eq. ([8]) that the behavior of the isovector incompressibility K\ should cor- 
relate closely with the NM symmetry energy S(p). We have plotted in Fig. the density 
dependence of Ki given by the numerical differentiation of the HF results for S(p) and one 
can see that the K\ value given by the (Asy-stiff) CDM3Yn interactions is positive over the 
density range p ^ p and gradually increases to 200 ~ 400 MeV at p approaching 0.6 fm ~ 3 
(the DU-onset density found with these interactions). On the contrary, the K\ value given 
by the (Asy-soft) M3Y-Pn and Gogny interactions is negative over the same density range 
and decreases linearly to about -1000 MeV at p approaching 0.6 fm ~ 3 . Since K(p) is actu- 
ally determined from the first derivative of NM pressure P{p) with respect to the density, a 
strongly negative isovector incompressibility K\{p) should indicate a decrease of P{p) in the 
transition from symmetric NM to the pure neutron matter. Indeed, one can see in Fig. [2] 
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FIG. 5: (Color online) Isovector part of the NM incompressibility given by HF calculations using 
different effective interactions. Encircled part shows K\ values near the saturation density po of 
symmetric NM, i.e., the K sym values given in Table HI See more details in text. 

that a decrease of P(p) in the pure neutron matter was found with the (Asy-soft) M3Y-Pn 
and Gogny interactions which pulls the calculated P{p) out of the empirical boundaries 



established by the HI collective flow data 



34|. Based on this discussion, we conclude that 



the behavior of isovector incompressibility K\(p) given by the (Asy-stiff) CDM3Yn interac- 
tions is more consistent with the HI flow data compared with K\(p) given by the (Asy-soft) 
M3Y-Pn and Gogny interactions. 

In the literature, the discussion on the isovector part of the NM incompressibility is very 
often made based on the K values estimated at the saturation density po of symmetric 
NM. It should be noted, however, that the saturation density of asymmetric NM decreases 
rather quickly with the increasing neutron-proton asymmetry 5, and the pure neutron matter 



(5 = 1) becomes unbound (see lower panel of Fig. [T] and also Fig. 2 of Ref. [17j). As a result, 
Po is no more a stable extremum in the NM energy curve and various expansions around 
it like (jSJ) might not be accurate for large neutron-proton asymmetries 5. For example, 
the second derivative of the approximated expression (j5J) for S(p) gives a purely parabolic 



18 



density dependence of the isovector incompressibility, K±(p) ~ K sym (p/po) 2 , which can 
strongly deviate from the exact HF result at high densities as shown in Fig. O In addition, 
the higher-order 0(S A ) term in Eq. fl5]) has been shown by Steiner 57J] to be important 
in determining the critical density for the DU process in neutron stars. In the vicinity of 
Po, the K\ values given by the HF calculation (encircled in Fig. ED are quite close to the 
corresponding K sym coefficients of the expansion (|SJ), i.e., K\{p — > po) ~ K sym . In the 



studies of the HI isospin diffusion 



391 ] or isospin dependence of the GMR excitation |60|, 



the asymmetry of the NM incompressibility around p was associated with the quantity 

K T = K sym — QL which has been confined by these data to K T —550 ± 100 MeV. The 

empirical K T value has been shown by a recent study of the neutron-skin thickness by 
I 

Centelles et al. [44j to be around K T « —500 ± 100 MeV. From Table [I] one can see that K T 
values obtained from the (Asy-stiff) DBHF, Vi ow k+CT, Hybrid and MDI (x=-l) results are 
in good agreement with the empirical values; the Asy-stiff CDM3Yn interactions give K T 
values of about -330 MeV which are somewhat above the upper limit of empirical data; the 
Asy-soft MDI (x=l), M3Y-Pn and Gogny interactions give K T values of about -400 MeV, 
right at the empirical upper limit. However, we note that a very recent mean-field study of 
asymmetric NM using a large number of Skyrme-type interactions by Chen et al. jsi] has 
found an optimum range K T m —370 ± 120 MeV and, thus, cast some doubt on previously 
adopted empirical K T values. In the context of our paper, this latest analysis seems to give 
preference to the (Asy-stiff) CDM3Yn interactions. Some further preference for the Asy-stiff 
interactions can also be made based on the empirical constraints on the J and L parameters 
established recently by the MSU group in the analysis of the isospin diffusion and ratios of 
neutron and proton spectra measured in HI collisions 33j. Namely, one can see in Tabled 
that only the L and J values given by the (Asy-stiff) CDM3Yn and SLy4 interactions are 
lying within the ranges of the double constraint deduced from the isospin diffusion data: 
L « 40 ~ 70 MeV and J « 30 ~ 34 MeV (see Fig. 1 of Ref. [33|). The (Asy-soft) Gogny 
and M3Y-Pn interactions give L values much lower than this limit and this indicates again 
that the (Asy-stiff) CDM3Yn interactions comply better with the HI data. We note here 
also a similar empirical range for the slope parameter L m 45 ~ 75 MeV established recently 
in a systematic study of the correlation between the neutron-skin thickness and symmetry 



energy 



62|. However, the neutron-skin thickness has been shown by Danielewicz 32| to be 



mainly sensitive to the surface part of the symmetry energy term in a more elaborate mass 
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formula for finite nuclei, while the extrapolation to the high-density behavior of the NM 
symmetry energy is based more on the volume term. Therefore, the empirical ranges for 
L, J and K T values deduced from the neutron-skin studies should be necessary reference 
points for any mean-field study of asymmetric NM but not sufficient constraints to restrict 
the behavior of the NM symmetry energy at high densities. 

Finally, we note that all the mean-field calculations discussed in the present work do not 
take into account the hyperon presence in the neutron star. The hyperon population has 
been estimated to make up about 18% of the neutron star matter and shown to significantly 



soften the EOS as well as reduce the limiting neutron star mass 63^ 



only the direct Urea process involving hyperons becomes possible 10J but also the proton 



64j. In this case, not 



fraction is significantly enhanced by the hyperon presence. For example, the proton fraction 
of a hyperon star having mass M ~ 1.5M is about 50% larger than that of a neutron- 
proton-lepton star of the same mass (see Fig. 5.28 of Ref. 64(). As a result, if we assume 
for simplicity a 50% rise in the proton fractions predicted, e.g., by the microscopic APR 
calculation or FSUgold model at p ~ 0.6 ~ 0.8 fm -3 (see lower panel of Fig. Hj), then the 
DU process is well allowed within these models. Concerning a typical Asy-soft interaction 
like the Gogny, M3Y-Pn or MDI (x=l), such a 50% increase of the proton fraction is still 
not enough to make the DU process possible. 



III. SUMMARY 



In the framework of the self-consistent HF mean field, we have studied the bulk nuclear 
matter properties predicted by two different sets (CDM3Y?i and M3Y-Pn) of the density- 
dependent M3Y interaction, SLy4 version of the Skyrme interaction as well as D1S and 
DIN versions of the Gogny interaction. The HF results for the NM symmetry energy and 
proton fraction in the /3-equilibrium are also compared with those given by the microscopic 
many-body studies (DBHF and APR calculations) using the bare NN interaction, and by 
the RMF studies using different parameter sets. 

We have concentrated our discussion on two main aspects: the NM binding energy and 
pressure in the symmetric NM and pure neutron matter, and the density dependence of 
the NM symmetry energy S(p) and the associated proton fraction (OQ). For the symmetric 
NM, the main conclusion is that all the effective NN interactions used here are more or 
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less consistent with the microscopic APR prediction and empirical pressure deduced from 
the HI collective flow data. For the pure neutron matter, the HF predictions for the NM 
binding energy and pressure show that the considered mean-field interactions are divided 
into two families which are associated with two different behaviors (Asy-soft and Asy-stiff ) 
of the NM symmetry energy at high densities, where only the Asy-stiff type interactions 
comply with the empirical NM pressure. These two families were shown to predict two 
different behaviors of the proton-to-neutron ratio in the /3-equilibrium which, in turn, imply 
two drastically different mechanisms for the neutron star cooling (with or without the direct 
Urea process). 

Although an ambiguity in the high-density behavior of the NM symmetry energy still 
remains due to the experimental evidences from HI studies favoring both the Asy-soft and 
Asy-stiff versions of the mean-field interaction, a comparison of the present HF results with 
the empirical constraints for the symmetry coefficient J and slope parameter L, given by the 
HI isospin diffusion data and ratios of neutron and proton spectra |33|] and the systematic 
study of the neutron-skin thickness 4J, I62J, seems to provide some evidence favoring the 
Asy-stiff type interactions. The Asy-stiff behavior is also predicted by the fully microscopic 



BHF or DBHF calculations 
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5n 



which include the hi ghe r-order many-body effects and 
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42]. 



three-body forces, and by the latest RMF studies 

A big puzzle remains why on the (nonrelativistic) mean-field level, a wide range of the nu- 
clear structure data can be consistently described only by using some Asy-soft type effective 
interaction like the famous Gogny forces or M3Y-Pn interaction. In each case, the chosen 
parameter set for the effective NN interaction depends strongly on the nuclear structure 
and/or reaction data under consideration 
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61j and there could be a plethora 



of systematic uncertainties in different choices for the mean-field interaction which are not 
under control and can lead, in particular, to the distinct soft and stiff scenarios. In any case, 
the ability of available structure and/or reaction data to constrain the EOS for high-density 
neutron rich NM on the mean-field level is still limited and the most interesting challenges 
are lying ahead. 
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